Exact reduced dynamics for a qubit in a precessing magnetic field and in the contact 
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The two-level quantum system (qubit) in a precessing magnetic field and in contact with a heat 
bath is investigated. The exact reduced dynamics for the qubit in question is obtained. We apply 
the approach based on the block operator matrices theory, in contrast with the standard methods 
provided by the theory of the open quantum systems. We also present the solution of the Riccati 
operator equation associated with the Hamiltonian of the system. Next, we study the adiabatic ap- 
proximation for the system under consideration using quantum fidelity as a way to measure validity 
of the adiabatic theory. We find that in the weak coupling domain the standard condition that 
ensures adiabatic behavior of the spin in the magnetic field also guarantees its adiabatic evolution 
in the open system variant of this model. Therefore, we provide the explicit example of the open 
quantum system that satisfies the adiabatic approximation firstly formulated for the closed quantum 
systems. 
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I. INTRODUCTION 

The most established and useful time-dependent two- 
level quantum system, is perhaps, the one that describes 
spin half particle (qubit Q) in a precessing magnetic 
field [1]. The Hamiltonian of this system in its basic 
variant is usually written in the following form 



ilQ(i) = -7S-B(f), S = *t/2, 



(1) 



where 7 is a constant, called the gyromagnetic ratio (the 
specific value of 7 is irrelevant in context of our discus- 
sion), (T = (cr^,CT*', cr^) and cr*, [i = x,y,z) are the stan- 
dard Pauli matrices. The magnetic field B(t) is assumed 
to have the form 



B(t) = Bi cos (ujt) ix + Bi sin (ujt) iy + BqCz 



(2) 



where ix, Sy, e^ are the unit vectors along x, y and z axis, 
respectively. Bq, Bi define the amplitudes of the mag- 
netic field and w is the frequency of the rotation. Using 
the equations above one can also rewrite the Hamilto- 
nian (1) in a more readable form: 



Hg^t, P) = paz + a [02 sin (wt) + cti cos (wi)) , (3) 

where for the sake of convenience the abbreviations a ~ 
iwi = —^7^1 and (3 = ^wq = —^jBo were introduced. 
In the case of no coupling with the external envi- 
ronment (heat-bath) the exact form of the evolution 
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operator UQ{t) and hence the density matrix pQ{t) = 
UQ{t)pQUQ{ty for the model (3), can be derived in an 
elegant and simple manner (see e.g. [I, 2]). This prob- 
lem is so common that it can be found in almost every 
modern textbook on quantum mechanics. However, if 
the aforementioned coupling is present, the exact form 
of the density matrix has not yet been derived. 

If one allows the system to interact (not necessarily to 
exchange the energy) with the environment then it be- 
comes so called open quantum system [■->] . Its time evolu- 
tion is not unitary any more because of the decoherence 
process [4, -5]. Nevertheless, it may always be described 
by trace preserving (TP) and completely positive (CP) 
quantum operation (or quantum channel, see e.g. [6]) 
Tt : pq ^ Pgit), Pq(0) = pg. Unfortunately, finding 
its exact form is almost impossible in most cases, espe- 
cially for the systems described by the time-dependent 
Hamiltonian. 

Naturally the question arises: how difficult this prob- 
lem actually is? Recently, the efforts to answer this ques- 
tion were made [7, 8] . It was shown that from the math- 
ematical perspective this task is at least as complicated 
as the problem of resolving the Riccati algebraic equa- 
tion associated with the Hamiltonian defining the model 
in question. Moreover, from the block operator matrices 
theory point of view obtaining the exact reduced dynam- 
ics for any system where decoherence is present (beyond 
dephasing phenomenon) [9, 10] is as complex as resolving 
this equation. Of course, the complexity of the analysis 
strongly depends on the particular choice of the environ- 
ment responsible for the decoherence process. 

First, we show that all the difficulties mentioned above, 
including the problem of resolving the Riccati equa- 
tion, can be overcome for the so-called spin environment. 
Therefore, we provide the exact reduced dynamics for the 
open system version of the model (1) in the case when 



the system is immersed within the spin-bath. Next, as an 
apphcation, we investigate the adiabatic approximation 
for the open system constructed in this way. Moreover, 
we show that in the weak couphng domain the condition 
ensuring the adiabatic evolution of the system (1) also 
leads to the adiabatic behavior of its generalization. 

The layout of this paper is as follows. In Sec. (II) 
we briefly review the concept of the block operator ma- 
trix and the procedure allowing one to diagonalize it. 
In Sec. (II A) we introduce the Riccati operator equa- 
tion and discuss its connection with the Hamiltonian 
specifying the model. In Sec. (Ill) the definition of the 
model is given and the exact reduced dynamics is derived. 
We also indicate possible applications of the model. In 
Sec. (V) we investigate the adiabatic approximation. Fi- 
nally, Sec. (VI) contains the final remarks and the sum- 
mary of the paper. 



II. BLOCK OPERATOR MATRICES 
APPROACH 

A. Partial trace and reduced dynamics 

Before we can go any further and present the main 
result of this paper let us first discuss the procedure al- 
lowing one to obtain the density matrix pqij:) for a two- 
dimensional system (qubit) using the block operator ma- 
trices approach. 

Let Hqe be the Hamiltonian of the closed (qubit and 
environment, Q+E) system. We will assume its following 
form 



Hqe = Hq®Ie+Iq®He + H„ 



(4) 



where Hq, He are the Hamiltonian of the qubit and 
the environment respectively, while Hint represents the 
interaction between the systems. Hamiltonian (4) acts 
on Hqe — "Hq (8) He space, where Hq and He are the 
Hilbert spaces for the qubit and the environment, respec- 
tively. The state PQ{t) of the qubit at any given time t 
can be computed using the following formula (for more 
details see e.g. [11] ) 



PQ{t)^TTE{UtPQ®PEUl) 

= Tt{pQ), 



(5) 



where Ut = cxjp{—iHQEt) is the evolution operator of 
the total Q + E system (we work with units fi, = 1). By 
TrB(-) we denoted the partial trace. In the literature the 
quantity pQ (t) is commonly called the reduced dynamics. 
In this paper we will be also referring to it as the solution 
of a given model. In Eq. (5) the assumption that pqe = 
PQ (^ Pe was made, i.e., no correlation between Q and 
E are present initially (see e.^., [12, 13] and references 
therein). 



Since Hq = C^ and the isomorphism C^ ® He = 
He® He holds, one can think of the operator Hqe (and 
of any given operator A acting on Hq (g) He) as a 2 x 2 
block operator matrix [14]. Let us assume that A acts 
on He ®He- Thus, we can write A = [Ay], where oper- 
ators Aij, {i,j = 1, 2) act on He- In this block operator 
matrices nomenclature the procedure of calculating the 
partial trace Tr^(-) can be written as 



TMA) 



TiAn TrAi2 
TrA2i TrA22 



(6) 



where Tr(-) denotes the trace on Hq. Eq. (6) intro- 
duces the concept of the partial trace in a very sim- 
ple and, what is the most important, intuitive way. It 
is also possible to define the partial trace locally, i.e., 
Tte{X(»Y) = XTt{Y), where as mentioned before Tr(-) 
stands for the trace operation on the space on which 
the operator Y acts on. From Eq. (6) one can learn 
that the partial trace is the operation which transforms 
operator matrices into the "ordinary" matrices. More- 
over, the partial trace is a linear operation going from 
B{Hqe) space to the i?(C^) space. Furthermore, from 
the block operator matrices perspective one can easily 
understand why dealing with the open systems is much 
more complex (from the mathematical viewpoint) then 
dealing with the closed systems. The reason is that the 
procedure allowing one to calculate the reduced dynamics 
is not as straightforward as Eq. (6) might indicate. In- 
deed, to compute the partial trace (5) one needs to write 
the evolution operator C/t in a 2 x 2 block matrix form. 
The last operation requires diagonalization of the block 
operator matrix Hqe, which is much more complicated 
than the diagonalization of the standard 2x2 matrix. It 
leads to the problem with solving the Riccati algebraic 
equation analyzed below. 



B. Riccati equation and block diagonalization 

With every Hermitian block operator matrix of the 
form 



R 



A B 
St C 



(7) 



where A, B, C are the Hermitian operators acting on 
Hei one can associate the following Riccati equation [15] 



XBX + XA-CX -B^ = 0, 



(8) 



where X G B{He) is the operator to be determined. 

If the solution X of the Eq. (8) exists, then it may be 
used to diagonalize operator matrix (7). The procedure 
that allows us to do so is as follows (for details see [8] 
and references therein) 



1 _\A + BX Oe 

where the similarity matrix Ux is given by 



(9) 



Ux 



X 






(10) 



We wish to emphasize that in most cases a general 
method of finding the solution for the Riccati equation 
does not exist. Therefore, resolving the Eq. (8) with re- 
spect to X is a very difficult task. Moreover, even if the 
solution is known, there is still a problem of computing 
the inversed operator matrix {Ux)~^- Formally, if the 
operator X is normal, i.e., [X, X^] = 0, then this matrix 
is given by the formula 



Hilbert space T-Le on which the Hamiltonian (13) acts 
is given by iV-fold tensor product of C^ spaces, i.e., 

N 

7i ^ (^ C^ . Therefore for any ji < N the operator a^ 

71=1 

is understood as tr^ = /2 (8) . . . (X) cr^ Cg) . . . (X) /2 , where a^ 
is the standard Pauli matrix acting on C^ and I2 is the 
2x2 identity matrix on that space. In our model the 
coupling of a qubit with the environment is governed by 
the Ising-typ Hamiltonian [IG, 17], namely 



N 



Hi, 



2_j9n<7ni 



(14) 



where .g„ are the coupling constants. We assume that the 
bath is initially in the Gibbs thermal state at a temper- 
ature T, i.e.. 



Ux^ = {Ie+XX^)-^ 



Ie X^ 
-X Ie 



(11) 



Otherwise, the form of the inverse operator U^ becomes 
more complicated. Fortunately, the Eq. (11) is sufficient 
for our analysis. Notice that it may be difficult to handle 
computationally the operator such as (I^ + XX'^)^^. In 
general the operator X does not need to be Hermitian or 
unitary. 

Perhaps those are the reasons why the approach to 
open systems based on the operator matrices theory has 
not attracted too much attention. In our opinion, how- 
ever, the operator matrix perspective can offers a better 
understanding of the open quantum systems and might 
give some new results in this area. 



III. EXACT REDUCED DYNAMICS 

A. Model 

We consider the model defined by the following time- 
dependent Hamiltonian 



HQEit, P) = Hgit, P)®Ie+Iq®He + H„ 



(12) 



where Hglt^P) is given by Eq. (3) and it represents a 
qubit in the magnetic field (2). The environment is com- 
posed of N independent and non-interacting spin 1/2 
particles. The Hamiltonian He of the bath is assumed 
to be of the form 



JV 



He = ^uJnK, 



(13) 



where a;„ and ct^, (n = 1,...,A'^) are the frequencies 
and the Pauli matrices for nth qubit, respectively. The 



P£ = Z-icxp(-i/B/A:T), 



(15) 



where Z = Tr(exp(— iJ^/fcT)) and k is the Boltzmann 
constant. 

The model described by the Hamiltonian (12) with He 
and V specified by the Eqs. (13)-(14) in the case of the 
a ~ (i.e., dephasing and static magnetic field case) was 
investigated both in the context of the approximation 
methods in the open quantum systems theory [Ki] and 
the capacities of the quantum channels [17]. Thus, for 
the detailed discussion and possible applications we refer 
the reader therein. 

The quantum system discussed here can by thought of 
as the generalization of the model introduced in the men- 
tioned papers to the case when the energy is exchanged 
between the systems and the magnetic field depends on 
time. It is worth mentioning physical problems that make 
use of this model. First of all, the Hamiltonian (12) may 
pose a useful prototype describing any spin "resonance" 
phenomenon, like for instance Nuclear Magnetic Reso- 
nance (NMR) [18, 19]. In such picture the spin 1/2 par- 
ticle is the open system, which evolution in time we wish 
to describe. The spin-bath models the influence of the 
other spins on the open system in question. Finally, the 
rotating magnetic held is used to induce the resonance 
and to control the system. 

The quantum devices of the future, like for instance 
quantum computers ['JO] will be probably built from the 
components that are composed of a large amount of the 
qubits (spins, in particular) e.g., quantum memory [21], 
quantum register [22], etc. The model we are investigat- 
ing may serves as a simple prototype allowing one to trace 
the evolution in time of a single qubit of the register or 
memory. The magnetic field may be then applied to pro- 
gram the device. Those are just a few potential applica- 
tions of the simple theoretical, time-dependent spin-spins 
model (12). 



B. Exact solution 

To derive an exact reduced dynamics of the system 
governed by the Hamihonian (12) we first simphfy this 
problem to a time-independent one. In order to accom- 
phsh this we use recently proved theorem that says if (for 
more details see (A) and also [S]) 

i) the total system is defined by the Hamiltonian (12), 

ii) the interaction term Hint between the systems 
takes the form f{a^) (g) V, where / is an analytical 
function of a^ and V is a Hermitian operator, 

iii) pt (/3) represents the solution of the system defined 
by the Hamiltonian Hqe{0, (3) (note that Hamilto- 
nian Hqe depends on /3, hence so is pi), 

then the reduced dynamics rjt for the model governed by 
the Hamiltonian (12) can be obtained using the following 
simple formula 



where Ux is given by the Eq. (10) and X is the solution 
of the Riccati Eq. (8), which in the present case takes the 
form 



aX^ + XH+ - H-X - alE = 0. 



(21) 



The quantity Hd represents a diagonal form of the oper- 
ator matrix Hqe- According to the Eq. (9) it reads 



Hd = 



H+ + aX Qe 

Qe H^ - aX^ 



(22) 



Note that for a = the solution (possibly not the only 
one) of the Eq. (21) is given by Xq ~ Qe- It is obvious, 
since in that case Hqe is already in the diagonal form, 
ie., Hqe = Hd ~ d\a,g{H ^ , H J) - In order to obtain the 
solution for a ^ more subtle investigation, provided in 
the next subsection, is needed. 



T^,=Vtpt{Peff)Vl /3e// := /? - ^/2, (16) 

where the unitary matrix Vt is given by (diag-diagonal) 



Vt 



diag I e 



-iujt/2 ibjtl'T 



(17) 



Equality (16) states that having the reduced dynamics pt 
one may easily obtain the solution we are interested in 
simply by introducing the effective parameter Peff, then 
replacing (3 by /3e// and finally performing the unitary 
transformation (17). 

The interaction term defined in the Eq. (14) satisfies 
the requirement of the theorem above, thus one should 
restrict the analysis to the Hamiltonian Hqe{P) = 
Hqe{0, P)- The later can be easily rewritten as a 2 x 2 
block operator matrix, namely 



HqeW) 



H+W) 
alE 



alE 



(18) 



where 



N 



^±(/3)-EK<±^^^)' P-Pi^^ (19) 



and LJ^^ :~ cj„ ± 5„, for each n < N- Henceforward, the 
explicit dependence of any quantity of parameter (3 will 
be omitted until the solution pQ{t) is obtained. 

The evolution operator for the total system generated 
by the Hamiltonian (18) reads 



Ut = cxp {-iHQEt) 

= Ux cxp {-iHd) Ux^ , 



(20) 



1- Solution of the Riccati equation for the a 7^ case 

First, one can observe that [H_,H^] = 0. Since the 
solution X of the Eq. (21) is a function of the operators 
H±, then it also must commute with those operators. 
The Riccati equation in question can be simplified to the 
more compact form, namely 



aX^ + 2VX -alE = 0, a^ 0, 



(23) 



where we introduced 



V 



N 

E 



9ncr„ 



Ph 



(24) 



Let us assume that the eigenvalue problem for the oper- 
ator V can be (easily) resolved. If this is the case we can 
write the solution X of the Eq. (23) in a manageable form 
using the spectral theorem for the Hermitian operators. 
By doing so we obtain that X = f{V), where 



/(A) 



VA2 +a^-X 



a 



, Xea{V). 



(25) 



Wc can also represent the solution f{V) in an equivalent 
way, i.e. 



f{v)= J2 fW\^)(M- 



(26) 



xe<j{v) 



Note that f = f* , i.e., / is a real function for any value 
of the parameters a and /3. Here, by criV) we denoted 
the spectrum of V . We want to emphasize that the rep- 
resentation (25) or (26) of the solution X of the Eq. (23) 



is useful computationally if the eigenvalues of V can be 
computed. We show below that this is the case. Note 
also that X =::XUSV = V^. 

Let i = iii2---iN be a binary (i„ = 0, 1, for n < N) 
expansion of an integer number i £ [0, 2^ — 1]. Clearly, 
the set of all eigenstates \i) = \ii) (g) \i2) (8 ... (8) Iin) forms 
an orthonormal basis in He- Moreover, because of the 
following equality CT^|i„) = (— l)*"|i„) we have 



N 



V\i) = J2{gr.{-iy-+^) 



(27) 



ES), 



0<i<2 



N 



1. 



Thus, the eigenvalue problem for V has been solved. 
Therefore, the solution of the Riccati Eq. (23) X reads 



N 



H±IO = EK(-i)'"±'^) 



(30) 



Et\i), 0<i<2 



N 



1. 



Observe also that E+ - E[ ^ 2Ei, for < i < 2^^ ~ 1, 
as one may expected. In the \i) basis the density matrix 
Pe can be expanded as 



Pe^J: (e-^'/'^^/z)\^){^\ 



i=0 



(31) 



where fii = S„'J-'n(~l)'" f-re the eigenvalues of He and 



X=J2 fiE^i^ 



(28) 



i=0 



At this point two remarks should by made. The first 
one is that the operator (28) depends on the parameter 
a, i.e., X = X{a) and X{a) — ;> Xq as a goes to 0. To 
see that this statement holds we rewrite the Eq. (25) in 
the form 



/(A) 



V^ 



Aecr(F). 



(29) 



1 



Clearly, /(A) —S- as a ^- 0, unless A = 0. The later 
means that the solution X of the Riccati Eq. (23) is a 
continuous function of the parameter a, including the 
a = value. From the physical point of view it means 
that one can control energy exchange between the sys- 
tems in a continuous way. Note that if a = then the 
energy transfer is absent. This scenario can be accom- 
plished by taking limits a — > with the final solution we 
are about to give. Note also that since for every a ^ 
and each A S <j{V) the eigenvalues /(A) of the operator X 
we constructed are positive. This solution is a positively 
defined operator on He. 

The second remark is that there exists, at least one 
more operator function that satisfies the Eq. (23), namely 
the one specified by the following formula: /2(A) ~ 
(— a/A^ + a^ — X)/a for A e <^{V). This solution rep- 
resents negatively defined operator. The choice between 
the different solution of the Eq. (23) is not arbitrary and 
it has significant influence on the further analysis. In- 
deed, if one decides to choose the second solution: /2(A) 
then one may meet a serious difficulties during the exam- 
ination of the dephasing phenomena, since /2(A) — > —00 
as a — > 0. 

At the end of this subsection let us note that one can 
resolve the eigenvalue problem for H± as easily as for the 
operator V . Indeed, we may readily verify that 



2. Total system evolution 

We are now in a position to give the explicit and man- 
ageable form of the evolution operator of the total system 
Q + E. First, let us note that the similarity operator ma- 
trix Ux can by written as 



Ux= Y. U^®\^){i\, 



(32) 



i=0 



where matrices Ui are given by 



t/z 



f 1 -.fiE.) 

[fiE^) 1 



(33) 



This follows immediately from the Eqs. (10) and (28). 
We also used the resolution of the identity I^; in the \i) 
basis. From the expansion (32) one can readily see that 
the inversed operator {Ux)~^ takes the form 



iUx)-'= E(^^)"'®N)(*i 



(34) 



where {Ui)~^ stands for the inverse of the matrix Ui. 
Since for every integer i wc have dct{Ui) = 1 + f'^{Ei) ^ 
0; thus the inverse exists. It is also easy to see that 
([/i)-i = dct{Ui)Ul and UiUj = l/dct{Ui). Because of 
the last equality, the operator (32) can become unitary 
by a simple rescaling procedure, i.e., Ui — >■ \/detUiUi. 
Furthermore, according to the Eq. (20) the diagonal form 
Uf := exp(— iiJrft) reads 



Uf 



E Uf{t)®\^){^\ 



(35) 



where abbreviations 



Ufit) = diag(e-*(^-^+"^'(^'»*,e-'(^'"-"/(^'»*), (36) 

were introduced. 

Finally, by combining all the results together we obtain 
the evolution operator Ut of the total system: 



2"-l 



Ut= J2 Urit)®\i){i 



(37) 



where 



form of the initial state Pqe- This assumption, however, 
is not essential and does not lead to the limitation of the 
analysis. The evolution operator (37) was written in a 
highly manageable form, and by that we mean it can be 
applied to any given initial state. Of course, if initial 
correlations arc present, then the Eqs. (39), and (40) 
no longer hold, but the reduced dynamics pqit) can still 
easily be obtained. However, we will not focus on this 
issue in the current paper. 

So far, we omitted the explicit dependency on /3 
in the introduced quantities. Note, however, that if 
Ui{t) ~ Ui{t,f3), then /9q(<) = pQ{t,(i) and in view of 
the Eqs. (16), (40), (17) and comments following them, 
we have 



um = u,ufit)ul 

_J_ ( etit) + er(t)fi Met{t) - er(t)) 
^ l + ff \.n (4 it) - 67 (0) e+ {t)fi + er (i) 

(38) 

and ef{t) := exp{^i{E^ ± af,)t) with /, = f{E,). Let 
us keep in mind that UiUj = Iq for each integer i. 



3. Reduced dynamics 

We begin with the assumption that pg and pE are ar- 
bitrary density operators for the open system and its en- 
vironment, respectively. Furthermore, we consider only 
factorable initial states pqe of the closed system Q + E, 
i.e., PQE = Pq ® Pe- Using this assumptions and the 
explicit form of the evolution operator Ut given in the 
Eq. (37), we have that the density operator PQE{t) ~ 
UiPqeUI at any given time reads 



m 



J2 M,,it)r,QM,jit)\ 



(41) 



where tjq [rjQ = pq) is an arbitrary initial state. The 
Kraus matrices Mij (t) are given by 



M.jit) = 5,j^VtU,{t,(i,}f). 



(42) 



The last two equations are the main result of the cur- 
rent paper. It is worth mentioning that the model with a 
similar properties was studied in [23], where the Authors 
obtained the operator sum representation, yet the Kraus 
operator introduced therein involved the time chronolog- 
ical operator. 

From the Eqs. (41) and (42) one can learn that the 
time evolution of the qubit interacting with the rotat- 
ing magnetic field (2) and in contact with the environ- 
ment (13) - (14) can be described using two quantum 
channels, i.e., PQ(t) — T^ o T^(pq). The first channel, 
that is T/ is given by 



PQEit) = J2 U^it)pQU,it)^ ® \i){^\pe\j){j\. (39) 

hj=0 

Tracing out the last equation over the environment de- 
grees of freedom, we obtain the reduced dynamics (note 
PqW = Tr£(pQB(i))) 



T/ := Y. mm^'hit)^ 



(43) 



The second one is defined by the simple unitary ope- 
ration, namely T^ = Vt(-)Vj. The matrices Mi are spec- 
ified by 



Pq(0 = Y. P^Mt)pQU^{t)\ 



(40) 



i=0 



where pi = {i\pE\i). By introducing matrices Kijit) 
such that Kijit) := Sij y/piUi{t) , the equation above can 
be written in a more familiar form, that is pgit) = 
J2ij Kij{t)pQKij{ty . Moreover, it is easily to verify that 
J2ij Kijit) Kij{ty = Iq. Thus, the Eq. (40) can be 
thought of as operator sum representation of the state 
Pqit). The matrices Kij{t) are the Kraus operators. The 
result (40) is general, i.e., it holds for the arbitrary states 
Pq and pE. Nevertheless, we made assumption about the 



Miit) = y/^,cxp{-iH^t), 



(44) 



where the Hamiltonian Hi is the generator of the unitary 
evolution Ui{t,f3eff) (sec (38)). It is easy to compute 
that 



„ _ (E^- cj/2 a 

^' ^ ' -E, + uj/2 



a 



(45) 



where Ei is the eigenvalues of the operator V and it was 
defined in the Eq. (27). The channel T^ may by thought 
of as the convex combination of the unitary channels 



(note Tt{lQ) = Iq). This type of channels are known 
as the random unitary channels [24] and they often ap- 
pears in the study of the pure decoherence. 

Surprisingly, the matrix above and therefore the Kraus 
matrices (44) does not depend on the function f{Ei). 
However, to compute matrices cxp(— i_ffit) that speci- 
fied the Kraus matrices (44) we need to diagonalize the 
Hamiltonians Hi. The function f{Ei) is hidden in the 
similarity matrix Ui given by the Eq. (33). Observe 
that the mentioned diagonalization procedure is based 
on the Riccati equation (matrices Ui are composed with 
the eigenvalues of the solution of the Riccati equation). 
This approach differs from the standard method based 
on the characteristic equation. This is a new kind of the 
diagonalization procedure, so called Riccati diagonaliza- 
tion [2.51. 



This is the well-known formula for the (unitary) evolu- 
tion operator of the system identified by the Hamiltonian 
HQ{t,f3) = Hqit). It is the solution of the differential 
equation idtU{t) = HQ{t)U{t). Note, that the result (47) 
corresponds to the Floquet decomposition of the unitary 
operator generated by the periodic Hamiltonian. Note 
also that the evolution operator (47) can by written in 
a more readable form, especially useful for studying the 
adiabatic approximation (see Sec. (V)), i.e. 



U{t) = VtCxp{-iHt), 



(48) 



where the matrix Vt is given by the Eq. (17) and H has 
the form 



IV. SPECIAL CASES 



H 



a 



a 



(49) 



Now, we look into some common cases that may arise 
during the examination of the model (12). We will show 
that in these particular situations the result (41) can be 
simplified to the well-know expressions. 



A. No coupling with the bath 

In the case when the qubit under consideration has no 
coupling with its environment, i.e., g„ = for all n < N 
we have 

a) Ei = f^eff, hence 

b) fi = f ^ i\z\ - I3eff)/a, where 

c) z := a + if3eff. 

Moreover, Ef = Qi ± (3, where Qi = X]n'^n(^l)'"- 
Using this simplification we obtain that Ui{t, ^eff) = 



-iUit 



U{t), where 



U{t) 



+ Pe 



2A\z\t 



-2ifsm{\z\t) 



-2i/sin(|z|f) /2e-*l^l* 
I2 cos{\z\t) — i 



„i\z\t 



^ J^j.AsHW) (46) 



~ I2 cos(|2|f) — ia ■ rism{\z\t) 
= exp(— i|z|(T • Ctt), 

with tl := n/||ri|| and Q. := (2/,0, 1-/^). As a result. 



M,,{t) 



-iQit 



— 5ijyfp'iVtU{t), where we dropped phase factor 
Therefore rjt = U {t)riU [t^ (as one may antici- 



pated), where 



U{t) = cxp(-iwtcr"'/2) exp(-i|z|(T • Q.t). (47) 



One can easily see that this matrix can be obtained di- 
rectly from the Eq. (45) if one set 5„ = for n < N . 



B. Dephasing 

If the constant a equals zero, then the Hamilto- 
nian HQE(t,f3) becomes time-independent. Moreover, 
in this circumstances [Hq, Hint] = hence the open 
system does not exchange energy with its environment, 
i.e., the pure decoherence occurs. It poses no prob- 
lems to check that C/i(t,/3e//) = e-*"'*Vj''' exp(-i£:icr^t). 
Therefore, the Kraus operators take the form Mij{t) — 
5ijy/piexp{—iEi(j^t). Just like before the phase factor 
g-iQit .^g^g drooped. In agreement with (41) the form 
above of the Kraus matrices leads to the following man- 
ifestation of the density matrix rjt 



Viiit) = ?7ii, Vi2{t) = Yl 



Pie 



i2_Eit 



?7l2 



(50) 



Naturally, 7]22{t) == 1 - r]ii{t) and r]2i{t) = V2iit)* ■ Note 
that the coherence C{t) :— |?7i2(i)| docs not decay expo- 
nentially (or anyhow for finite TV) in the long time regime, 
but manifests oscillating behavior. 

From the considerations above one can easily learn 
that there is a substantial difference between the situ- 
ation when a transfer of the energy from the qubit to its 
bath is not present and the case when the qubit is not 
coupled to the environment. Let us also observe that in 
the model we study it is not possible to construct a sit- 
uation when the energy exchange in not present {a = 0) 
and yet the magnetic field is rotating. This follows from 
the fact that if a = then the Hamiltonian HQE{t,(3) 
does not depend on u. 
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C. Low temperature regime 

From the mathematical point of view the zero (low) 
temperature limit is the most subtle special case. This 
comes from the fact that most results obtained by using 
the known approximation methods in the open quantum 
systems theory, e.g., the Markovian or the singular cou- 
pling limit, cannot by extrapolated to this regime. Luck- 
ily, the exactly solvable models are unhampered by this 
difhculties. 

To derive the exact reduced dynamics in the low tem- 
perature domain we begin with rewriting the density op- 
erator pe in a more suitable form in contrary with the 
one given by Eq. (31). According to the result obtained 
in [IG] one may write 



N 



PE = 



(Iq + (3„a') , 



(51) 



where /3„ = tanh(— w„//cT). From equality (51) one may 
readily see when T — >■ 0, then 



;5^n(l + /3"(-ir") 



6io, 



(52) 



since /3„ — > 1 as T ^ 0. All this means that in the low 
temperature domain the heat bath is its ground state, 
i.e., Pe ~ \0e){0e\- As an immediate result we obtain 
the following form of the reduced dynamics (41) 



lim rjt 



Wt-qWi 



(53) 



where Wt := VtUo{t, (3^ff). Interestingly, the evolution of 
the open system in question is unitary, which is rather un- 
usual. This is a direct consequence of the fact that for ev- 
ery integer i the operator Ui{t) specified by the Eq. (38) is 
unitary. This means that in the low temperature regime 
the dissipative properties of the heat bath (13) - (14) are 
frozen. 



D. Equal Couplings and Frequencies 

Henceforward and without essential loss of generality 
we assume that a;„ = f7 and gn = 9 for n < N, i.e., 
couplings constants and the frequencies of the spins of 
the bath are equal. In this situation Ei = g{N — 2k) + /3, 
where k is the Hamming weight of the integer number i 
(i.e., a number of nonzero element in a binary expansion 
of i). Since there are ( j,) integer number i G [0,2^ — 1] 
with the same Hamming weight k, the channel (43) takes 
the form 



with Mk{t) given by Eq. (44). Note, in this case pk (see 
Eq. (51)) simplified to the form 



Pk^^{i+5r-Hi-5f 



(55) 



where 5 := ian^-Vl / kT) . 



V. APPLICATION: ADIABATIC 
APPROXIMATION 

The adiabatic theorem [26, 27] for the closed quan- 
tum system specified by the Hamiltonian H{t) in its ba- 
sics variant {i.e., discrete and no degenerate spectrum 
a{H{t))) states that if H{t) varies slowly (for rigorous 
meaning of that see e.g., [2n] and references therein) and 
if the system is initially prepared in one of the eigen- 
statcs of H{0), say |V'n(0)), then at any given time t the 
probability \{tpn{t)\'4'{t))\'^ of finding it in the eigenstate 
\ipn{t)) (an eigenvector of H{t)) is equal to 1. Therefore, 
one may easily introduce the quantity F{t) that measures 
the validity of the adiabatic approximation, namely 



F{t) = Tripit)p.^,{t)), 



(56) 



where p^{t) = \^pt){M and pit) = UtpiO)UJ with p{0) = 
|V'o)(V'o|- Here, l-^t) is the eigenvector of H{t) for i > 
and Ut represents the (unitary) evolution operator for 
the closed system. If the system is evolving adiabati- 
cally then F{t) = 1. Note that definition of the quantity 
F{t) coincides with the standard definition of the quan- 
tum fidelity F{p{t),p^{t)) between pure states p{t) and 
p-^,{t) [2!)]. Since this is a two dimensional case it also 
coincides with the definition of the super fidelity [M), 31]. 
Observe that the Eq. (56) is very useful since it can be 
easily applied both to the closed and to the open quan- 
tum systems. In the last case we need to replace the 
equality p{t) = Utp{0)Uj representing the evolution of 
the closed system by the channel describing evolution in 
time of the open system. More precisely, p{t) = Tt{p{0)), 
where Tt is TP-CP quantum operation. In the current 
section we use the Eq. (56) in order to investigate behav- 
ior of the model (12) in the adiabatic regime. To make 
the paper self-contained we briefly discuss the adiabatic 
approximation for the system (1) first. 



A. The closed spin system case 

For the sake of simplicity we put /3 = f3o cos ip and 
a = /3o sin ip for certain (p and /Sq. This assumptions lead 
to the following form of the Hamiltonian (3) 



N 



Tl=E 



fc=0 



N 



Mk{t){-)Mk{t)^ 



(54) 



HQ{t)-Po(^,T/^^ 



e *"* sin If 
— cos ip 



(57) 



It can be easily verified that the matrix above has the 
eigenvalues E^ — ±/3o- The corresponding eigenvectors 
l^f ) are given by 



\^t) 



cos((^/2) 
e*"* sin(^/2) 



(58) 



The density matrix p+(i) = \il}^){tpf\ takes the form 



P+{t) 



cof?{Lpl2) ie-''"*sin¥J 
ie''^*sin(y9 sin^((^/2) 



(59) 



Observe that p+{t) = Vtp+{0)V^ , where Vt is speci- 
fied in the Eq. (17). According to the Eq. (48) the 
density matrix p{t) at any given time t reads p{t) — 
VtU{t)p+{Q){VtU{t))\ thus the fidelity (56) takes the 
form 



F(f) = Tr(C7(i)p+(0);7(i)V+(0)) 



= \\U{t)pm\\ 



(60) 



_Fi 



where || • \\f is the Frobenius (Euclidean) norm, i.e., 
||A|||,:=Tr(AAt). 

It is well-known that for the system governed by the 
Hamiltonian (57) condition that guarantees the adiabatic 
evolution is /3o 3> w. This statement is very intuitive, it 
means that the magnetic field rotates slowly, in compar- 
ison with the phase of the state vector (for more detail 
discussion see e.g., [32]). We can also see this from the 
Eq. (60). Indeed, if one introduces the adiabatic param- 
eter X := oj/2f3o, then in agreement with (48) and (60) 
we have 



Fit) = 1 - 



1 + 2:' 



■sm'^{n{x)t), 



(61) 



where n{x) :— /^oVT+j?. Without loss of generality 
we set if = tt/2 in the equation above. The parameter 
X measures how slowly the magnetic field rotates in /3o 
units, thus in the adiabatic limits (i.e., x — > 0) the sec- 
ond term in the Eq. (61) can be neglected and therefore 
F{t) c. 1. 



B. The open system case 

By analogy to the closed spin case discussed above one 
may write the fidelity (56) for the open system (12) as 



Fit) = Tr{T^oTl{p^{0))p4t)) 



N 

E 

A:=0 



N 



PkFkit), 



(62) 



where Fk{t) := \\Uk{t)p^{0)\\%. The channels T/ and 
T^ are specified in the Eq. (43) and comments bellow. 



respectively. The unitary matrix Uk takes the form Uk ~ 
exp(— iTJfet), with Hk given by the Eq. (45). If we put 
(fi = -K 12 as before, then 



Fk(t) = 1 - 



^l 



^ + 4 



sm'^{n{xk)t), 



(63) 



where Xk := G{N - 2k) /Po - x. Using this and Eq. (62) 
we obtain that F{t) = 1 — R{t), where R{t) reads 



N 



R{t) = Y. 



k=0 



,PkSm^{n{xk)t). (64) 



From this equation one can readily see that in the adi- 
abatic domain {x — > 0) we have R{t) ^ and thus 
F{t) < 1. Therefore, the standard condition that leads to 
the adiabatic behavior of the closed system (1) does not 
guarantee that in the case of the open system model (12) 
it holds true as well. However, if one additionally assumes 
that coupling with the environment is weak, in compar- 
ison with the energy split between the states |0) and |1) 
of the qubit in question, that is to say if G/(3o ^ 1, then 
ccfc < 1 (for finite N) and R{t) ~ 0, thus F{t) ~ 1. 



VI. SUMMARY 

In this paper we investigated a qubit in contact with 
the spin environment and interacting with a rotation 
magnetic field. The considered model was constructed 
under a set of assumptions which allow for its exact treat- 
ment. We hope that despite the mathematical character 
the results of the paper may serve as starting point for 
further investigations. The exact models can not only 
provide reasonable approximate description of real sys- 
tems, as it is the case for the pure dephasing, but are 
often used as a basis and inspiration for various improve- 
ments [■-'>?>]. Although the paper is mainly focused on 
mathematical aspects it also includes the example of a 
natural application of the model in question. This exam- 
ple relates our model to important problems of physics 
such as e.g. the problem of the adiabatic approximation 
for the open quantum systems [34] or the adiabatic quan- 
tum computation [35]. 

We provided the exact reduced dynamics for the sys- 
tem mentioned above. In contrast to the standard meth- 
ods available in the theory of the open quantum systems 
our approach was based on the block operator matrices 
theory. In particular, we resolved the algebraic Riccati 
equation associated with the Hamiltonian defining the 
model under consideration. We wish to stress out that 
the method we used in the current paper, although ap- 
plied to the particular model, is general. Nevertheless, 
its usefulness relies on the ability of solving the Riccati 
equation. At the present time it is a very difficult task, 
even for the simple systems. 
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At this point one can ask: how relevant is it to assume 
a spin environment, instead of a bosonic one, to solve 
the problem exactly? Is it possible to do so for bosonic 
bath? One may pose a more general question. What 
is a connection between the possibility of obtaining the 
exact reduced dynamics of the qubit in question and the 
Hamiltonian specifying the bath? This problem was al- 
ready addressed in [s] . The results of this paper as well 
as the analysis we carried out suggest that this problem 
is at least as difficult as resolving the Riccati equation 
associated with the total Hamiltonian. 

Furthermore, we studied the adiabatic approximation 
for the model in question. It was shown that the standard 
condition that guarantees the adiabatic evolution in the 
the case of the closed systems is not valid for the open 
system generalization. This is not an unexpected result. 
It is interesting, however, that the aforementioned con- 
dition does ensures the adiabatic behavior of the open 
system under consideration in the weak coupling limit. 
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Appendix A: Proof of the Eq. (16) 

In order to prove the equation (16) let us note that the 
Hamiltonian (12) satisfies the following condition 



H{t,P)=e^'^'H{P)e 



-iKt 



(Al) 



where K — — ^0-3 (g) Ig. This can be easily proven us- 
ing the Baker-Campbell-Hausdorff formula [36]. As was 
shown, in [2] every quantum system with Hamiltonian 
H{t,(3) satisfying (Al) for some Hermitian operator K 



evolves 

Ut{l3) = e»^*e-*^=/^(^)*, i?e//(/3) := H{f3) + K. (A2) 

Note that in general [H{I3),K] ^ and therefore 
[Heff{f3),K] 7^ 0. In our case, from equation (12) we 
learn that H{f3) = {(icr^ + acri) ® Ie, thus 



HeffiP) == (/3ct3 + aCTi) ® Ib - -(73 ® Ie 



((^ 



0-3 + aai 



(A3) 



From equations (A2) and (A3) we have 



um 



iKt 



C/t(/3-f), 



(A4) 



where Ut{l3) is the evolution operator generated by 
H(t,fi). Let pt{P) and fjt be a density operator for the 
closed system Q + E associated with the Hamiltonian 
iJ(/9) and H{t,f3) respectively in some arbitrary time t. 
Let us also assume that po{(3) = fjo = p. Using equa- 
tion (A4) one can easily see that 



f,t = Ut{f3)pUHf3) 



(A5) 



'Oe 



-iKt 



where we introduced Vt = e''^*. 
will show that ii Ai, A2 <E BCH i 
operator matrix of the form Ai ~ 
B = [B,j] e B{n ® n) then 



To end the proof we 
D •H) are a 2 X 2 block 

Ai ®Ie, {i = 1,2) and 



Tve{AiBA2) = AiTve{B)A2. 



(A6) 



Equation (A6) follows from the linearity of the trace Tr(-) 
operation and the definition (6) of partial trace. Note 
that Vt = Vt®lEi where Vt is given by the equation (17), 
thus taking partial trace of the equation (A5) and us- 
ing (A6) we obtain (16) with Vt given by (17). 
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